Band offsets and stability of BeTe/ZnSe (100) hetero junctions 
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We present ab-initio studies of band offsets, formation energy, and stability of (100) hetero- 
junctions between (Zn,Be)(Se,Te) zincblende compounds, and in particular of the lattice-matched 
BeTe/ZnSe interface. Equal band offsets are found at Be/Se and Zn/Te abrupt interfaces, as well 
as at mixed interfaces, in agreement with the established understanding of band offsets at isovalent 
heterojunctions. Thermodynamical arguments suggest that islands of non-nominal composition may 
form at the interface, causing offset variations over ~ 0.8 eV depending on growth conditions. Our 
findings reconcile recent experiments on BeTe/ZnSe with the accepted theoretical description. 

73.40.Lq, 73.61.Ga. 
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On the basis of experiments as well as theory, it is 
commonly accepted that at isovalent semiconductor in- 
terfaces the band offset is almost independent of the 
local atomic arrangement, except iiju the presence of 
heterovalent interlayers or antisites.Elo This result was 
originally established for the common-ion systems, and 
later generalized to the whole class of isovalent het- 
erojunctions. Investigations on no-common-ion systems 
such as InP/GalnAs and InAs/GaSb confirmed that the 
band offset is independent on the atomic-scale inter- 
face arrangement, despite -the different interface com- 
position and local strain.|jQ Remarkably, these findings 
have found a rationale wiiiin the linear response the- 
ory (LRT) of band offsets,E[l3 which also aoeounts for the 
composition-dependent local strain effects.Ll 

Only at heterovalent junctions does the band offset 
depend crucially on interface morphology, due to the dif- 
ferent chemical valence of the atomic constituents, and it 
is fully explained within LRT. So far, the maximum vari- 
ation experimentally detected amojunts to 0.6 eV, and it 
was observed at ZnSe/GaAs (100).l3 

Controversial findings have been reported for the iso- 
valent lattijae-matched BeTe/ZnSe (100) interfaces. In a 
first paperjj a valence band offset (VBO) of 0.9 eV was 
deduced from the luminescence spectca of BeTe/ZnSe 
quantum wells. In a subsequent work,E3 thin BeTe films 
grown on ZnSe (100) were investigated by XPS: unex- 
pectedly, two widely different values, 0.46 eV and 1.26 
eV, were measured in different growth conditions, and 
interpreted as due to Se- and Zn-terminated substrates. 

Should this interpretation be confirmed, this result 
would be a) the first case of morphology-dependent band 
offset at isovalent interfaces; b) the largest VBO variation 
(0.8 eV) ever observed at semiconductor heterojunctions; 
c) a clear violation of the LRT of offsets. 

To help the interpretation of the experimental data, 
here we investigate the band alignment and the thermo- 
dynamical stability of BeTe/ZnSe(100), and other related 
junctions among (Zn,Be)(Se,Te) compounds, using first- 
principles density-functional-theory calculations. These 
have proven to be a highly reliable tool in predicting off- 



sets at semiconductor-semiconductor interfaces. 

Pseudopotential plane-waves caleulations are per- 
formed using ptiie VASP code, til with ultrasoft 
pseudopotentialsc3 for Be and Zn (including Zn 3d states 
in the valence). iTie generalized gradient approxima- 
tion (GGA-PW9lll3) to the exchang e-correlation func- 
tional is used. All relevant properties of the binary 
(Zn,Be)(Se,Te) compounds are converged at a cutoff of 
23 Ry. As shown in Table ||, the experimental lattice 
parameters are reproduced within 1%. The abrupt in- 
terfaces were modeled in periodic boundary conditions 
by 24-atom slab supercells. For mixed or reconstructed 
interfaces, 48-atom supercells with a total of 24 atomic 
layers and 2 atoms per layer were used. Brillouin zone 
integration was performed on a 6x6x2 Monkhorst-Pack 
mesh. The in-plane (substrate) lattice parameter in the 
supercell calculation was chosen to be asub=5.697 A, i.e. 
the average of the theoretical bulk lattice parameters of 
BeTe and ZnSe. 

For each interface, the supercell structure was fully 
optimized, as it is mandatory to obtain realistic results, 
since Zn-Te (Be-Se) bond lenghts differ by about ± 10% 
from those of bulk BeTe and ZnSe. We relax ionic po- 
sitions and cell parameters until forces below 0.05 eV/A 
and stress along the (100) direction lower than 0.5 Kbar 
are obtained. The VBO is then computed following the 
approach described in Ref. |[ The comparison of the 
VBO calculated for the ideal unrelaxed and for the op- 
timized structures confirm the importance of structural 
optimization, whose effects on VBO amount to about 0.4 
eV. 

Two kinds of abrupt interfaces are possible at 
BeTe/ZnSe(100) heterojunctions: the Zn/Te interface, 
characterized by the sequence of atomic planes . . . - 
Be-Te-Bc-Te-Zn-Se-Zn-Se- . . . , and the complementary 
Be/Se interface, with the stacking sequence ...-Te-Be- 
Te-Be-Se-Zn-Se-Zn- .... In supercells with periodic 
boundary conditions, interfaces are always present in 
pairs, and they may be chosen (for the present orien- 
tation) to be different or identical, depending on the 
atomic filling of the supercell. We consider supercells 
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both with identical, symmetry-equivalent interfaces (the 
first two columns of Table |l|) and asymmetric supercells 
with different interfaces (the third column of the Table, 
marked " asym" ) . The comparison of the results allows to 
reduce the numerical uncertainty in the estimate of the 
VBO. Supercells with symmetry-equivalent interfaces ex- 
hibit a non-ideal c/a ratio due to local strain in the in- 
terface regions, whereas in supercells with different and 
complementary interfaces the local positive and negative 
strains nearly compensate, and c/a is close to the ideal 
unstrained value. The average VBO for the abrupt re- 
laxed (100) interfaces is 0.52 eV, and all values fall within 
a range of ~ 30 meV, which is of the same order of magni- 
tude of the numerical uncertainty of the calculations 
10 meV), in analogVjto the results for other no-common 
ion heterojunctions.B^O The equivalence of the two inter- 
faces occurs also for the ideal unrelaxed cases. Therefore 
LRT is valid also for this system, despite the large chemi- 
cal differences between the constituting compounds (e.g., 
BeTe and ZnSe bulks have a very diflferent ionicity, 0.34 
and 0.59 on the Garcia-Cohen scalgiil). The validity of 
LRT is also confirmed by the fact that, for the unrelaxed 
case, the VBO for the abrupt (110) interface is close to 
the average of the two different (100) terminations. The 
VBO's between differently oriented, relaxed interfaces 
shows larger differences; but this is not in contrast with 
the LRT picture, because local interface strains depend 
not only on composition but also on orientation. From 
all the above results, we definitely rule out the possibility 
that the difference of about 800 meV between the VBO's 
measured at BeTe/ZnSe (100) interfaces can be pmply 
ascribed to chemically different abrupt interfaces. lij 

As a further check, we studied some prototypical (100) 
non-abrupt interfaces, restricting to either anion- or 
cation-intermixed cases (antisites are generally energeti- 
cally unfavorable in II- VI compounds). In particular we 
consider c(2x2) reconstructed interfaces with one mixed 
layer of either Be and Zn atoms (fourth column in Ta- 
ble H) or Se and Te atoms (fifth column in the Table). 
Again, the VBO is independent of the interface local 
atomic arrangement within 10 meV. Therefore, we con- 
clude that the VBO at BeTe/ZnSe (100) heterojunction 
does not depend on the interface local atomic arrange- 
ment, thus confirming previous evidence for isovalent in- 
terfaces, and the general predictions of LRT. 

The DFT-GGA VBO values are not directly compa- 
rable with experimental data, since spin-orbit coupling 
and self-energy effects on bulk bands are not taken into 
account in this type of electronic structure calculations. 
The spin orbit splitting is 0.96 eV for BeTe and 0.40 eV 
for ZnSc; including a posteriori the ensuing correction to 
the VBO reported in Table ||, we obtain an estimate of 
about 0.7 eV. Many-body corrections to valence band top 
edges, still excluded from this estimate, are not available 
for these compounds to our knowledge. Typical values 
for these corrections are of order 0.1-0.2 eV,E£l so that a 
final theoretical estimate could be close to the experimen- 
tal value of 0.9 eV reported in Ref. ^ and, incidentally, 



to the average of the two values of Ref. ^ However, we 
stress that the corrections to the DFT-GGA VBO val- 
ues are bulk quantities, and thus they affect the absolute 
value of the VBO, but not at all the relative compari- 
son among the values for different cases considered here. 
Therefore the main result of our calculations, i.e. the in- 
dependence of the VBO on interface composition, is fully 
valid. 

According to our calculations, the VBO of 0.9 eV re- 
ported in Ref. ^ could correspond to several possible in- 
terface compositions, including either Zn/Te or Be/Se 
abrupt terminations or mixed interfaces. However, some 
suggestions about the actual structure of the interface 
comes from a thermodynamic investigation of interface 
stability. We find that in thermodynamic equilibrium 
abrupt interfaces of either kind are favored over the in- 
termixed ones. We define the interface formation energy 
per unit of sectional area in the most general case as 
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where E"^l^ is the total energy of the supercell describ- 
ing the interface, and the /i's and iV's are the chemical 
potentials and number of atoms of the various elements 
involved. At equilibrium the chemical potentials of the 
elements and total energies of the condensed phases are 
related by 
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The formation energy of abrupt interfaces is easily seen 
to be a function of the difference between Zn (or Se) and 
Be (or Te) chemical potentials. Indeed, using Eqs.|l[ the 
formation energy for abrupt Zn/Te and Be/Se interfaces 
reads 

O TTiZn/To _ pZn/Tc ,^ ,,BoTc t\t ,,ZnSe / ,Zn ,,Bc\ 

r, j-iBc/Sc _ pBc/So ,r ,,BcTc at ,,ZnSc, /,,Zn ,,BoN 
^%orm - ^tot -AtcM " Asc M + (^J. j, 



respectively. The range of variation of /i^" — /i^'^ is 

,,Zn ,,Bc ^ ,,Zn-bulk ,, Bo-bulk a tjBoTo 
/i — MSM ^ Llll , 

,,Zn ,,Bc \ ,,Zn — bulk Be — bulk i a rjZnSc 

where AH-^ is the formation entalpy for compound 
X. Mixed-interface supercells are instead stoichiometric 
(Asc=Azn=AznSc, and Atc=Abc=AboTo), therefore the 
formation energy is independent of the chemical poten- 
tials. The previous expression becomes 
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where now A's and /I's are referred to the bulk formula 
unit. The results, summarized in Figure |^, show that 
the Zn/Te abrupt interface is favored in high (^^" — M^°) 
conditions, and conversely the Be/Se abrupt interface is 
favored in low (/i^'^ — fjP°) conditions. Most interest- 
ingly we find that, unlike the case of heterovalent junc- 
tions, the present isovalent abrupt interfaces are always 
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favored over the mixed ones for the whole range of admis- 
sible chemical potentials. This behavior was alreadyj^re- 
dicted for the III-V isovalent GalnP/GaAs interface|13 so 
we suggest that this preference for abrupt interfaces may 
be generally valid for any isovalent heterojunction. 

Thermodynamics further gives key indications (at least 
as far as equilibrium energetics is concerned) on the pos- 
sible origin of different offsets measured in real samples 
in particular growth conditions. Islands of a priori un- 
expected composition, such as BeSe or ZnTe, may form 
during the deposition of BeTe on ZnSe: specifically one 
expects BeSe islands in Be-rich and Se-rich growth con- 
ditions, and ZnTe islands in Zn-rich and Te-rich condi- 
tions. In terms of band offsets, the idea is that these 
"hetero-islands" may in fact be the material effectively 
interfaced to ZnSe, and therefore largely determine the 
observed band offset. 

The idea of islands formation is suggested by previous 
experience with dopant incorporation in semiconductors. 
It was shown theoreticallytZl that rising the chemical po- 
tential of the Li acceptor up to its bulk value, Li incorpo- 
ration in ZnSe is preempted by the formation of a Li2Se 
surface phase— Indeed, heavy Li doping of ZnSe layers in 
MBE growthlia results in the formation of Li2Se islands 
on the ZnSe surface. In the present case the scenario 
is slightly more complex, as four chemical potentials are 
involved. We choose as reference the cation chemical po- 
tentials, both for convenience and because the cations 
are the mobile species; the phase diagram of the four- 
component interface system will thus be drawn in the 
{/i^", /i^'^} plane. The reactions leading to the forma- 
tion of an epitaxial compound on ZnSe at the expenses 
of BeTe are as follows: for ZnTe on ZnSe, 



Zn + BeTe ^ Be ZnTe, 
and for BeSe on ZnSe 

Be + ZnSe =^ Zn + BeSe. 



(2) 
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These reactions will occur exothermically if the reaction 
energy AE is negative; the latter energy is given for re- 
actions H and H by 



^gZnTc 
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respectively. In these relations, /i^^ is the total energy 
of bulk XY in the pseudomorphically strained geometry 
on ZnSc, as it results from the optimized XY/ZnSe inter- 
face supercell. Using these equations and the calculated 
values of the chemical potentials and compounds forma- 
tion energies, we determine the regions in the {/z^",^^''} 
plane where BeTe and ZnSe are unstable with respect to 
trasformation into ZnTe and BeSe. The phase diagram is 
represented in Fig|^. In region A, AE^''^° is negative and 
/\]T;ZnTc jg positive: therefore the formation of epitaxial 
BeSe through reaction |^ is energetically favored. In re- 
gion B, AE^°^° is positive and AE^""""" negative, hence 



epitaxial ZnTe is energetically favored over BeTe. In re- 
gion C, both the AE's are negative, hence both BeTe and 
ZnSe are unstable respect to decomposition into BeSe 
and ZnTe. According to this picture, at thermodynami- 
cal equilibrium BcTc/ZnSc interfaces are never stable and 
the following interfaces may locally form instead: refer- 
ring to Fig. |, BeSe/ZnSe in region A, ZnTe/ZnSe in 
region B, and ZnTe/BeSe in the (very small) region C. 
Our present result indicates that interfaces established 
in real BeTe/ZnSe samples might be locally closer to 
ZnTe/ZnSe in Zn-rich conditions and BeSe/ZnSe in Be- 
rich conditions, than to the nominal BeTe/ZnSe compo- 
sition. A direct consequence of this result which should 
be observable in experiment is the preferential formation 
of BeSe or ZnTe islands on ZnSe during the early stages 
of growth of a nominally BeTe-ZnSe interface. Our anal- 
ysis does not include growth kinetics effects, which may 
cause the (unstable) nominally-BeTe/ZnSe interface to 
actually form for chemical potentials in region C of Fig. 
^, where the thermodynamic driving force towards equi- 
librium (i.e. instability of BeTe/ZnSe) is smallest. 

We now discuss the key piece of information we are 
looking for, namely the VBO values for the various pos- 
sible interfaces. In calculating them, we use the same 
in-plane lattice parameter asub as in all previous calcu- 
lations (the substrate is unchanged), and carefully ac- 
count for bulk and interfacial strain effects. As in the 
BeTe/ZnSe case, the calculated VBO values are affected 
by a substantial absolute uncertainty due to many-body 
and spin-orbit splitting effects, here combined with split- 
tings coming from epitaxial strain. However, the rela- 
tive uncertainty in comparing the values for the different 
systems is much smaller, due to a partial cancellation 
of systematic corrections to the bulk band edges. The 
results are depicted schematically (also including the no- 
common-ion lattice-matched interface) in Fig. ||. Two 
points are relevant in the Figure: (a) the transitivity 
ruleEl holds within the numerical uncertainty of the cal- 
culations, confirming once again the validity of LRT for 
these systems; (b) the values of the VBO for the differ- 
ent systems differ at most by about 0.8 eV, the minimum 
value corresponding to the BeSe/ZnSe interface (Be-rich 
conditions) and the maximum to ZnTe/ZnSe (Zn-rich 
conditions). 

It is interesting to note that the maximum calculated 
VBO difference of 0.8 eV is the same as the one mea- 



sured in the two different samples in Ref. 10; in addition. 



in that experiment the maximum value was observed in 
Zn-rich conditions, and the minimum value in Se-rich 
conditions, in agreement with our findings. This match- 
ing suggests a possible correspondence hntwecn the lower 
(higher) experimentally measured VBOllj and the forma- 
tion of a BeSe/ZnSe (ZnTe/ZnSe) interface, although the 
absolute values of the calculated offsets are about 0.3 eV 
lower than the measured ones, because of the discussed 
unaccuracy of the theoretical estimate. 

In conclusion, we presented band offset calculations 
for a series of zincblende (100) interfaces between var- 
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ious (Zn,Be)(Se,Te) II- VI compounds. We also set up 
a thermodynamical phase diagram bearing on the sta- 
bihty of the various possible interfaces. Based on our 
results, we discussed recent experiments on BeTe/ZnSe 
interfaces, which showed a marked offset variation with 
growth conditions. Our conclusions are that: (i) the 
attribution of the two widely different measured VBO 
values to abrupt Se-terminated and Zn-terminated in- 
terfaces of the nominal BeTe/ZnSe heterojunction, as 
proposed in Ref. |l^, is incorrect, as well as any other 
attribution to mixed (reconstructed) interfaces, which 
have a composition-independent VBO; (m) conversely, 
strained interfaces between other (Zn,Be)(Se,Te) com- 
pounds shows a VBO which may differ up to 0.8 eV in 
the case of BeSe/ZnSe and ZnTe/ZnSe; (in) thermody- 
namics indicates that such interfaces may actually locally 
form in the deposition of BeSe on ZnTe and viceversa; 
(iv) interfaces between (Zn,Be)(Se,Te) compounds follow 
closely the linear reponsc theory predictions just as III- 
V-based systems. Although the problem require further 
investigation for a definite explanation, the comparison of 
experimental and our theoretical findings could suggest 
that observed interfaces may locally be not the nominal 
BeTe/ZnSe, but rather interfaces such as BeSe/ZnSe or 
ZnTe/ZnSe depending on the chosen growth conditions. 
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FIG. 1. Formation energy (in eV/(lxl) cell) of the abrupt 
and c(2x2)-reconstructed BeTe/ZnSe (100) interfaces as func- 
tion of the difference of the Zn and Be chemical potentials. 
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FIG. 2. Phase diagram for the (Be,Se,Zn,Te) common-ion 
interfaces. The stable interfaces are BeSe/ZnSe in region A, 
ZnTe/ZnSe in region B, BeSe/ZnTe in region C. 
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FIG. 3. Schematic of the valence band offsets for the vari- 
ous interfaces investigated. 



TABLE I. Theoretical and experimentaltj lattice param- 
eters (in A) of the zincblende compounds considered in this 
work. 

Compound BeTe BeSe ZnSe ZnTe 

Exp. 5.626 5.139 5.668 6.104 

GGA-PW91 5.661 5.174 5.734 6.182 



TABLE II. Valence-band offset (in eV) at different 
BeTe/ZnSe interfaces. The valence band top edge is higher 
in BeTe than in ZnSe. 
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